Long-range coherent coupling in a quantum dot array 

Floris R. Braakman 1 *, Pierre Barthelemy 1 , Christian Reichl 2 , 
Werner Wegscheider 2 , and Lieven M. K. Vandersypen 1 * 
1: Kavli Institute of Nanoscience, TU Delft, 
2600 GA Delft, The Netherlands and 
2: Solid State Physics Laboratory, ETH Zurich, 8093 Zurich, Switzerland 

(Dated: March 6, 2013) 

Abstract 

Controlling long-range quantum correlations is central to quantum computation and simulation. 
In quantum dot arrays, experiments so far rely on nearest-neighbour couplings only, and inducing 
long-range correlations requires sequential local operations. Here we show that two distant sites can 
be tunnel coupled directly. The coupling is mediated by virtual occupation of an intermediate site, 
with a strength that is controlled via the energy detuning of this site. It permits a single charge to 
oscillate coherently between the outer sites of a triple dot array without passing through the middle, 
as demonstrated through the observation of Landau-Zener-Stiickelberg interference. The long- 
range coupling significantly improves the prospects of fault-tolerant quantum computation using 
quantum dot arrays and opens up new avenues for performing quantum simulations in nanoscale 
devices. 



Nanofabricated circuits of quantum dot arrays provide an excellent platform 
for quantum information processing using single charges or spins [TrH]. In such 
a dot array architecture, the tunnel coupling between neighbouring dots plays 
an essential role. It governs the motion of charges between the dots, permitting 
derealization |5J over the dots and coherent oscillations between them J3J [6], 
and the same tunnel coupling is at the core of exchange-based quantum gates 
on spin qubits [H, El E] • 

Tunnel coupling falls off exponentially with distance, and all experiments on 
quantum dot arrays so far rely on nearest-neighbour couplings only. In addition, 
quantum dot arrays are typically constructed from one-dimensional segments 
since realizing two-dimensional arrays is challenging. These restrictions severely 
constrain the range of experiments possible in this system at present. Instead 
of having to repeatedly swap neighboring qubits down the chain, a long-range 
coupling would enable quantum gates between distant qubits in one step, thereby 
giving access to many of the benefits of a two-dimensional lattice. This would 
strongly reduce the requirements for fault-tolerant quantum computing [9l 110] 
and permit quantum simulation of phenomena that are otherwise inaccessible 
in this system, for instance involving frustration [llj or superexchange [12, 13J. 

The most common approach to realizing an effective long-range coupling is to use a 
quantum bus, as demonstrated for trapped ions [14J and superconducting qubits [To] . 
For quantum dots, such a bus has been proposed in the form of optical cavities [17] and 
microwave stripline resonators [T8H20] . For the latter, the first steps have been taken re- 
cently [2T1 122] . Furthermore, charge transfer through a channel connecting two distant 
quantum dots has recently been realized using surface acoustic waves that push electrons 
forward [231 123] - In this approach charge coherence is lost but spin coherence is expected to 
be preserved. 

As an alternative for creating long-range coupling of quantum dots, which does not require 
separate elements such as cavities or channels, a quantum bus has been proposed in the form 
of the continuum of conduction or valence band states [25] . Through a second-order process 
known as cotunneling, virtual occupation of these states can induce an effective coupling 
between distant quantum dots. Inspired by this scheme, we propose to create such long- 
range coupling by virtual occupation of discrete states of quantum dots located in between. 
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In this case only discrete levels participate in the cotunneling process. This permits a fully 
coherent process, in contrast to all existing measurements of cotunneling in quantum dots 
in which quantum coherence is quickly lost in the reservoirs (see [25] for a review). 

Here we demonstrate the coherent transfer of single electron charges between the outer 
sites of a linear array consisting of three quantum dots, in a regime where sequential trans- 
port through the middle dot is suppressed energetically. Using real-time charge detection 
techniques, we study the dependence of the rate with which electrons hop between the outer 
dots on the detuning of the middle dot levels. We observe a non-monotonous dependence 
that is characteristic of a coupling mechanism mediated by virtual occupation of the middle 
dot levels. We also control quantum coherent dynamics between the outer dots in the form of 
Landau-Zener-Stuckelberg (LZS) interference, induced by a process we dub photon-assisted 
cotunneling (PACT). 

A scanning electron micrograph of a device identical to the one used is shown in Figure 
la. Gate electrodes fabricated on the surface of a GaAs/AlGaAs heterostructure (see Sup- 
plementary Information) are biased with appropriate voltages to selectively deplete regions 
of the two-dimensional electron gas (2DEG) below and define the linear array of three quan- 
tum dots. In the array, only adjacent dots are connected through tunnel barriers. The left 
and right dots are also tunnel coupled to the left and right reservoirs respectively. Above 
the blue-shaded gate a charge sensing dot (SQD) is created, the conductance of which is 
sensitive to the number of charges on each dot in the array through capacitive coupling. For 
maximum sensitivity, the SQD is operated on the flank of a Coulomb peak. Furthermore, 
one of the SQD contacts is connected via a bias-tee to an LC-circuit so that the SQC con- 
ductance can be measured both by RF reflectometry [27] and in DC (see Supplementary 
Information). 

We operate the device in the few-electron regime: by sweeping the voltages on gates LP, 
MP and RP, the number of electrons on each of the dots of the triple dot array can be changed 
one by one. A time-averaged measurement of the differential DC-conductance through the 
SQD as a function of Vlp and Vrp maps out a cross-section of the three-dimensional charge 
stability diagram of the triple dot (Figure lb). The occupancy is denoted here as (n,m,p), 
corresponding to the number of electrons on the left, middle and right dot respectively. 

The strength of the four tunnel couplings can be tuned individually with the voltages 
on the B-gates. The tunnel rates between the outer dots and their respective leads are set 
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Figure 1: a, SEM image of a sample identical to the one used for these measurements. Dotted 
circles indicate quantum dots, squares indicate Fermi reservoirs in the 2DEG, which are contacted 
through Ohmic contacts. Both the current through (white arrow) and the reflectance of the SQD 
is monitored and used to determine the occupancies of the triple quantum dot. b, Numerical 
derivative (along Vlp axis) of the current through the SQD as a function of the voltages on gates 
LP and RP, mapping out a charge stability diagram of the triple dot in the few-electron regime. 
The (0,0,0)-(0,l,0) charging transition appears fragmented because of low tunneling rates from 
the reservoirs to the center dot. c, Real-time traces of the sensing dot reflectometry signal, taken 
at points L, R and C in Fig. lb, as indicated. We use a 50 kHz low-pass filter (Avens Signal 
Equipment AP220) to filter the reflectometry signal in order to obtain sufficient signal-to-noise. 

to the 100 Hz range. The tunnel rates between neighbouring sites are tuned to be much 
higher. The upper two panels of Fig. lc show real-time traces of the charge detector RF 
reflectometry signal taken at points R and L in the charge stability diagram. The first trace 
shows a single step, corresponding to the transfer of one electron from the left dot to the 
left reservoir, i.e. going from (1, 1,0) to (0, 1,0). In the second trace, three single-electron 
tunnel events are seen, once from the right dot to the right reservoir and twice the other 
way, i.e. alternating between (0, 1, 1) and (0, 1,0). The step size is larger than in the first 
trace, because of the closer proximity of the SQD to the right dot than to the left dot. 
Remarkably, when we go to point C in the charge stability diagram, the real-time trace 



(lower panel) not only shows steps corresponding to the slow tunneling between outer dots 
and leads, but also exhibits smaller steps that occur at a rate which is an order of magnitude 
higher. Since point C is at the boundary of the (1,1,0) or (0,1,1) regions in the charge 
stability diagram, the fast steps appear to correspond to single electron transfers between 
the outer two dots. This is consistent with the step size as well as with the observation that 
the mean times that the measured conductance is high or low are equal for zero detuning 
between (1, 1,0) and (0, 1, 1) (point C). Upon increasing or decreasing the detuning, these 
times quickly become unequal (see supporing online material). 

This tunneling between the left and right dots is at first sight unexpected, since in these 
measurements the center dot levels are far detuned from resonance with the outer dot levels, 
excluding sequential tunneling via the center dot (there is no charging line of the center 
dot present nearby in the charge stability diagram). Here we argue that these tunneling 
events are transfers of single electrons between the outer dots, via virtual occupation of the 
middle dot. This implies that electrons are transferred between the outer parts of the array, 
essentially without passing through the dot in between. 

This tunneling between remote dots can be seen for different charge configurations of the 
triple dot array. Here we focus on transitions between (1,1,0) and (0, 1, 1). In this case, two 
virtual pathways exist for the transfer: either a single electron moves first virtually from 
the left to the middle dot and then from the middle to the right dot, or an electron moves 
first from middle to right and then another electron moves from left to middle (note that for 
other charge configurations the situation can be different, for instance for tunneling between 
(1,0,0) and (0,0, 1), only the first pathway exists). As will be shown below, the existence 
of two virtual pathways makes the transfer rate depend non-monotonously on the detuning 
between the intermediate virtual states and the initial and final states. The dependence on 
the middle dot detuning is a key signature of this process. 
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Figure 2: a, b, Schematic representations of the cotunneling process in terms of the relevant 
electrochemical potentials in the linear dot array. The two panels illustrate the two possible 
pathways for cotunneling between 1 1 10) and 1 01 1) , as explained in the main text. 

The charge transfer is depicted schematically in Figures 2a and 2b. Since it involves 
removing an electron from one dot and adding it to another dot, we need to compare 
electrochemical potentials for the transitions between initial, intermediate virtual, and final 
charge states. Only ground-state to ground-state transitions are considered. We denote the 
various electrochemical potentials as: 

/ii(n, m,p) = E(n,m,p) — E(n — l,m,p) 

HM{n, m,p) = E(n, m,p) — E(n,m — l,p) 

Unin, m, p) = E(n, m, p) — E(n, m, p — 1) 
for the left, middle and right dot respectively. Here E(n,m,p) stands for the ground-state 
energy of a given charge configuration (n,m,p). For the transition back and forth between 
the states 1 110) and 1 1 1 ) to occur spontaneously, we require /xl(1,1,0) = /ijj(0, 1, 1) (we 
denote the lowest energy state with occupation (n,m,p) as \nmp)). Furthermore, in order 
to quantify the detuning between virtual states, and initial and final states, we introduce 
two parameters 8\ and 82- For charge transfer via virtual occupation of |020), the relevant 
detuning parameter is 81 = /iA/(0,2,0) — //l(1, 1,0) (see Fig. 2a). For charge transfer via 
1 101) , the relevant detuning is 82 = /zl(1,0, 1) — /za/(1,1,0) (Fig. 2b). Note that 8± and 82 
are related and cannot be changed independently, as Vmp increases 8\ by the same amount 
it decreases 82 ■ The total tunnel rate T is the sum of the tunnel rates via the two respective 
paths. It can be expressed as (see Supplementary Information): 

h \ 81 + 81 ) [l) 
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Here tn, t r i, U 2 and t r 2 are the tunnel coupling elements between 1 1 10) and |020), 1 020) 
and 1 01 1) , ]| 101 ) and |011), and 1 1 10) and 1 101) , respectively. The charge dephasing time 
T 2 is assumed to be much smaller than 1/r (T 2 is typically about 1 ns [H E] and therefore 
interference effects between the two pathways can be neglected). Equation 1 is valid as long 
as the four tunnel couplings and the detuning between 1 110) and |011) are small compared 
to 5i and 62- 

We experimentally verify the non-monotonous dependence of T on detuning 5i (and hence 
82) by stepping the voltage on gate MP, Vmp, and measuring the rate of tunneling between 
1 1 10) and 1 011) . Figure 3a presents three traces, each for a different value of Vmp- For 
the top trace the value of Vmp corresponds to small 5i and large 62, therefore the transfer 
proceeds mainly as depicted in Fig. 2a. For the middle trace, Vmp is set such that both Si 
and 62 are relatively large, resulting in a reduced, but non-zero tunnel rate, in agreement 
with Eq.l. Finally, for the lower trace, 62 is small and Si is large. In that case, tunneling 
proceeds mainly via the virtual process shown in Fig. 2b and the tunnel rate is higher again. 

For a quantitative analysis, we extract the tunnel rate T from real-time traces such as in 
Fig. 3a using the relation [28J T -1 = /(l — /) ((tl) + (tr)), with / the Fermi probability 
distribution. The times (tl) and (tr) are the average times between tunnel events spent 
in the left and right dot, respectively. We perform a threshold analysis of the real-time 
traces to obtain the distributions of tl and Tr. The value of / can be established using the 
relation / = (r L ) / ((r L ) + (t r )). Figure 3b shows measured values of T thus determined, for 
different values of the detuning, parametrized by 5\. The non-monotonous dependence is 
striking and is fit well by Eq.l (red curve), implying that the transfer indeed proceeds via 
virtual occupation of intermediate states. 

The hopping of electrons between the outer sites in the array indicates that an effective 
tunnel coupling is present between the left and right dot, which we call cotunnel coupling. 
We can express the strength of this cotunnel coupling as (see Supplementary Information) 

_ tnt r i tl2t r 2 /„\ 

We note that t co need not be small: for typical experimental values of the nearest-neighbour 
tunnel couplings of order 10 /zeV and for detunings of ~100 /ieV, we obtain t co ~ 1 /xeV. 
The cotunnel coupling enters the Hamiltonian in much the same way as the direct tunnel 
coupling between neighbouring sites and therefore many phenomena arising from direct 
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Figure 3: a, Real-time traces of the SQD reflectometry signal, taken at zero detuning between 
the outer dot levels, for three values of Vmp, corresponding to three values of 5\ and 82- b, Plot 
of the measured cotunneling rate T versus detuning 6%. The non-monotonous dependence is a 
clear indication that the transfer proceeds via cotunneling. This is corroborated by the fact that 
the measured data points can be fitted well with the predicted expression for T (red curve). To 
make the fit, we rewrite Eq.l as T = a/5f + b/(c — d5\) 2 , with a, 6, c and d positive constants. For 
the detuning axis, gate voltages are converted to energies using microwave-induced sidebands as 
an energy reference (see Supplementary Information). The error bars on the obtained values for T 
include errors associated with the threshold analysis of the real-time traces (low-frequency noise 
modulates the baseline signal, so the precise value of the threshold slightly affects the statistics) 
and sampling errors due to the finite number of transfer events per trace [28J (we sample over 100 
ms traces). Note that the use of a low-pass filter results in an overall underestimation of T. 

tunnel coupling have their counterpart in cotunnel coupling between remote sites. For 
instance, upon application of microwave excitation to a gate, direct tunnel coupling can 
give rise to photon-assisted tunneling (PAT). In complete analogy we can expect a photon- 
assisted version of the cotunnelling process described above, which we term photon-assisted 
cotunneling (PACT). 

Photon-assisted tunneling (PAT) is a well-described phenomenon in quantum dots and 
has been observed many times in single and tunnel coupled double quantum dots [HI 129] . 
Tunneling transitions of electrons between two detuned neighbouring dots can be made 
resonant by applying microwaves of a frequency matching the detuning, Sq = nhv (see 
Figure 4a). Here Eq is the detuning, v the microwave frequency and n an integer, showing 
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Figure 4: a, Schematic view of photon-assisted tunneling processes between different pairs of 
dots. Charges can be transferred from one dot to another when the detuning between the 
corresponding electrochemical potentials matches the photon energy. The left and middle panels 
correspond to PAT, the right panel corresponds to PACT. Note that similar resonances to the 
ones shown exist for negative detunings. b, Charge stability diagram in the same configuration as 
in Fig. lb, but now with microwave excitation (15 GHz) applied via a bias-tee to gate LP. The 
microwaves were chopped at the reference frequency of a lock-in amplifier and combined with a 
small amplitude modulation of the same reference frequency. The colorscale data is the numerical 
derivative (along Vlp axis) of the SQD signal acquired via the lock-in amplifier. Multiple 
sidebands develop where PAT or PACT occurs. 

that also multiphoton resonances are possible. 

We apply microwave excitation at v — 15 GHz to gate LP. The microwaves are chopped 
at the reference frequency of a lock- in amplifier (see Supplementary Information). The 
excitation introduces a number of sidebands in the charge stability diagram. Two sets of 
sidebands are due to conventional PAT between neighbouring tunnel coupled quantum dots. 
They are seen in Figure 4b near the | LOO) to 1 0L0) transition (point L) and 1 00L) to 1 0L0) 
transition (point R). The slope of these lines is such that eo = nhv is maintained. 

Near the | L LO) to |011) transition (point C), a different set of resonances develop. Based 
on their slope and location, we identify these transitions to occur via PACT, where sin- 
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gle electrons tunnel between the outer dots, now assisted by the microwave excitation. As 
expected for photon-assisted processes, these resonances appear at a detuning linearly de- 
pendent on the microwave frequency (see Supplementary Information). 

In order to get a strong PACT response, the tunnel couplings between neighbouring dots 
are set to a much higher value than in the real-time experiment descibed before. We can 
extract the nearest-neighbour tunnel couplings from the spacing between the PAT resonances 
along the detuning cLXIS clS db function of frequency j5j [29] (see Supplementary Information) 
and find 8.1 ± 0.4 /xeV for the tunnel coupling between right and center dot, and 12.3 ± 0.3 
fieV between left and center dot. 

Importantly, photon-assisted cotunneling allows us to demonstrate the coherent dynamics 
driven by the cotunnel coupling, through the observation of Landau-Zener-Stiickelberg (LZS) 
interference [30] - Coherent quantum dynamics in the form of LZS interference has been 
observed in a wide range of two-level quantum systems, such as electronic states of atoms 
and molecules J3U, E2] or superconducting devices [33] and spin states in double quantum 
dots [31] • In LZS interferometry, a two- level system is swept through an anticrossing of its 
levels at such a rate that a superposition of its ground and excited states is reached (see 
Figure 5a). Between two passings through the anticrossing (at times t\ and tz), the two 
parts of the superposition acquire a relative phase due to their difference in energy e, 



At the second passing through the anticrossing, the two paths in phase space will interfere. 
Destructive interference in the occupation probability of the excited state occurs for A0 12 = 
(2n + 1)tt, where n is an integer. 

In our experiment, the applied microwaves modulate the detuning e between /xl(1, 1,0) 
and fin(0, 1, 1). This takes the system back and forth through the anticrossing created by the 
cotunnel coupling (Fig. 5a). The amplitude of the microwaves is used to control the value 
of A0i2 [33]. Figure 5b shows a measurement of the SQD signal as a function of detuning 
Sq versus microwave power. Contrast against the background indicates a finite population 
in the excited state (the ground state configuration for positive and negative detuning is as 
indicated in the figure). We see that for both positive and negative detunings, the excited 
state population exhibits interference fringes when sweeping the power (moving along the 
vertical dashed line). The oscillations indicate that a coherent superposition of 1 110) and 




(3) 
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Figure 5: a, Schematic energy level diagram as function of detuning between 1 1 10} and 1 01 1) , 
displaying an avoided crossing due to cotunnel coupling. The red arrows represent the response of 
the system to microwaves modulating the detuning. Multiple passings of the avoided crossing 
results in quantum interference of the two paths, b, Numerical derivative (along detuning axis) of 
the lock-in signal of Isqd (as in Figure 4b) as a function of detuning and microwave power. 
LZS-interference fringes are clearly visible along both axes. 

|011) is created and maintained between subsequent passings through the anticrossing. 

Since the microwaves drive the system through the anticrossing periodically, also multiple 
crossings need to be considered. Consecutive cycles interfere constructively when the total 
phase difference accumulated over one complete microwave period AO tot equals 27m. As 
A@tot — hj e(t)dt = 27C6Q/hu, this causes peaks to appear for detunings £o,n = nhu. This 
can be clearly seen in Figure 5b, where along the horizontal direction ten fringes can be 
discerned. 

In summary, we have demonstrated an effective coherent cotunnel coupling between the 
outer dots in a triple quantum dot, which is mediated by virtual occupation of levels on the 
dot in between. The coupling strength can be controlled through the detuning between the 
relevant middle and outer dot levels and agrees well with theoretical predictions. 

The long-range tunnel coupling may be used as well for realizing spin exchange gates 
at a distance in one step. When the intermediate site is itself occupied by an electron, 
its spin affects the strength of the cotunnel coupling due to the Pauli exclusion principle, 
introducing correlations between the middle spin and the outer spins (i.e. it realizes a three- 
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qubit gate). First schemes for avoiding these correlations have been worked out, enabling 
direct long-range spin exchange also in this case [35J. Furthermore, the cotunnel coupling 
can be extended to include multiple intermediate sites. Long-range coupling thus provides 
a new approach for operating quantum circuits based on quantum dot qubits, which eases 
the requirements for fault-tolerance. The cotunnel coupling observed here also gives access 
to a new range of phenomema with interacting spins, such as superexchange [T2l [T3] and 
frustration [11], which can serve as a starting point for quantum simulations. 

Methods 

The experiment was performed on a GaAsZAlo.25Gao.75As heterostructure grown by 
molecular beam epitaxy, with a 85 nm deep two-dimensional electron gas with electron 
density of 2.0 x 1CT 11 cm -2 and mobility of 5.6 x 1CT 6 cm 2 V _1 s _1 at 4K. The metallic (Ti- 
Au) surface gates were fabricated using electron beam lithography. The device was cooled 
inside an Oxford AST Sorb dilution refrigerator to a base temperature of ~55mK. In order 
to reduce charge noise |36j, the sample was cooled down while applying a positive voltage 
bias on all gates, ranging between 200 and 350 mV. The magnetic field, as well as the bias 
across the linear triple quantum dot were set to zero throughout the experiment. Gates 
LP and RP were connected to homebuilt bias-tees, enabling application of DC as well as 
high-frequency voltage bias to these gates. RF reflectometry of the SQD was performed 
using an LC circuit matching a carrier wave of frequency 193.35MHz. The power of the 
carrier wave arriving at the sample was about -84 dBm. The reflected signal was amplified 
using a cryogenic Quinstar QCA-U-219-33H amplifier and subsequently demodulated using 
homebuilt electronics. A Stanford Research Systems SR830 lock-in amplifier was utilized in 
some of the measurements. In these measurements, a square wave modulation of amplitude 
2mV, before 16 dB attenuation, was applied to gate LP at the lock-in reference frequency 
of 3412Hz. For the microwave measurements, this square wave was combined with the mi- 
crowaves, chopped at the same reference frequency. The microwaves were generated by an 
Agilent E8267D microwave source. 
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Supplementary Information for 
Long-range coherent coupling in a quantum dot array 



I. ADDITIONAL CHARGE STABILITY DIAGRAMS 
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Figure SI: a, Numerical derivative (along horizontal axis) of the conductance of the SQD as a 
function of Vlp and Vpp, measured in the opaque barrier regime of Fig.lc. b, Larger range 
measurement, displaying the numerical derivative (along horizontal axis) of the SQD conductance 
as a function of Vlp and Vpp, this time in the more open barrier regime of Figs. lb and 4b. 



The real-time traces of Figure lc and Figure 3a were taken in a different configuration 
of gate voltages than where the charge stability diagrams of Figure lb and the PACT data 
of Figure 4b and Figure 5b were taken. Figure SI. a shows a charge stability diagram in 
this first regime, which was tuned such that all barriers (between dots and between dots 
and reservoirs) were quite opaque. Note that charging lines of the middle dot are not very 
visible, since it is charged at a very low rate due to the presence of multiple high barriers 
between this dot and the reservoirs. Figure Sl.b shows a charge stability diagram in the 
more open barrier regime of Fig.4b, this time for a larger range of gate voltages than shown 
in the main text. 
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II. REAL-TIME TRACES FOR DIFFERENT DETUNINGS BETWEEN OUTER 
DOTS 
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Figure SI: Real-time traces of the SQD signal for different detunings between (1, 1,0) and 
(0, 1, 1), taken for <5i ~300//eV. From top to bottom the detunings are: 49^eV, 12^eV, -24/ieV 
and -61/xeV. 



Supporting evidence that the tunneling events shown in Fig.lc are between the outer 
dots is given by the traces of Figure S2. The different traces are taken for different values of 
the detuning e between (1, 1, 0) and (0, 1, 1). From top to bottom, the detuning is increased 
from negative to positive values. In the upper trace, mainly (0,1,1) is occupied. In the 
center traces, charges go back and forth most rapidly between the outer dots and in the 
lower trace, mostly (1, 1, 0) is occupied. 



17 



III. COTUNNELING: EFFECTIVE TUNNEL COUPLING 



H = 



(4) 



The Hamiltonian describing the experiments can be expressed in the basis = 
|110), \4 0) ) = |011), |4 0) ) = |020), |4 0) ) = |101) as: 

f-e/2 t n t l2 \ 
t r \ t r2 

til t r i Si 

\ tn t r2 S 2 J 

We perform a unitary transformation of this Hamiltonian, to express it in the eigenbasis 
of its first-order perturbation: H' = WHU, where, as e « Si >2 : 

( 1 tn/Si t l2 /S 2 \ 

o i t rl /s 1 t r2 /s 2 

-tn/Si -t rl /S 1 1 
\-ti 2 /S 2 -t r2 /S 2 1 / 

In this new basis, it is sufficient to only consider processes between the two first states 



U 



(5) 



(i)\ 



1 110> - ^|020) - ^|101) 

Si s 2 



|^) = |011) -^|020) -^|101) 
Si s 2 



(6) 
(7) 



Processes involving mixing between ijjp or ip^ and the other states are of third order in 



U/Sj and can therefore be neglected. We can then reduce the Hamiltonian to: 
H = 



hi 
Si 



h'2 tqtrl _|_ tl'jtrl 
62 Si 82 



Ultrl I tl2tr2 e \r\ 

Si " r 8 2 2 Si 




This Hamiltonian then simply expresses an effective tunnel coupling t c 



+ 



(8) 



tl2~tr2 
82 



between states whose detuning has been renormalized. In the presence of microwaves, we 
can write the detuning as e — > Eq + Aexp(iuut), and the Hamiltonian becomes: 



H = 



(9) 



—s f /2 — Aexp(iut) t co 

t co e' /2 + Aexpiiut) 

which is exactly the Landau-Zener-Stuckelberg Hamiltonian, so all the physics of the LZS- 
interference can be directly applied to describe photon-assisted cotunneling in the triple dot 
array. 
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IV. LOWER BOUND ESTIMATION FOR CHARGE T 2 



The width w of the PACT resonances can be used to establish a lower bound on the 
charge dephasing time T 2 , according to the relation 



'-^w+^h (10) 

With the measured value w = 4.47 GHz, we find a lower boundary for T 2 of 224 ps. From 
PAT measurements between the left and middle and right and middle dots we extract lower 
bounds for T 2 of 394 ps and 338 ps respectively. 

V. CALCULATION OF THE REAL-TIME TRANSITION RATE 

The relatively slow charge transition process in Figure 3 can be understood by means 
of a density matrix description. The evolution of the system is described by the quantum 
Liouville equation, -£ = —i/h[p } H], where H is the Hamiltonian and p the density matrix. 
Given an effective tunnel coupling t co between two resonant states, the Hamiltonian is: 



H 



t co 

trn 



Adding phenomenologically decoherence occuring on a timescale T 2 , we get the system 
of equations: 



d 

dt Pll = 


it 

h 


(P21 - 


P12) 


d 

dt p22 = 


it CO 

IT 


(P21 ~ 


P12) 


d 

dt Pl2 = 


Pl2 

T 2 




(p22 ~ Pll) 


d 

dt p21 = 


P21 

T 2 


it 


(P22 - Pll) 



We do not include relaxation channels, since processes inducing charge transitions other 
than the cotunneling are much slower than both cotunneling and decoherence. Introducing 
N = p 22 — pn and P = p 2 i — pi 2 , the equations reduce to 



dt h 

d „ P 2it co Ar 

— P = h — -N 

dt T 2 n 
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We are interested in the rate with which charge moves back and forth between the outer 
dots. In our system, the initial conditions describe a pure state, where the electron is in 
the left dot: N(0) = —1, P(0) = 0. The rate at which charge tunnels to the right dot is 

dp-njdt = 7,dN/dt, where 



For a coherent process (t co » 1/T 2 ) we obtain that a charge oscillates between the two 
dots at a rate linearly dependent on the coupling t co : 



In the time traces of Fig. lc and 3a, the barriers are tuned so that tunneling is extremely 
slow (t co « I/T2), and the transition rate becomes 



As the timescale of the measurement is much longer than the decoherence time, the charge 
transition rate is therefore given by: 



In the experiments described in the manuscript, two paths ( 1 1 10) — y 1 101) — y |011) and 
1 1 10) — y |020) — y 1 01 1) ) contribute to the coupling t co . For the slow, real time measurements 
of the tunneling rate, these two contributions add up incoherently and we have 



where ti = y/tnt r i (t 2 = y/t^t^) describes the coupling through 1 101) (|020)), and 5± and 
5 2 are the respective detunings. 

VI. CALIBRATION OF THE DETUNING BETWEEN MIDDLE AND OUTER 
DOT LEVELS IN FIG. 3B 

We wish to determine the detunings 5\ and 82 between the intermediate and initial/final 
electrochemical potentials for any set of gate voltages. To do so, we first convert gate 





dt 



2t 2 co T 2 (l - e^ 2 ) 
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Figure S3: a, Numerical derivative (along horizontal axis) of the SQD lock-in signal (as in 
Fig. 4b) as a function of Vlp and Vrp, while microwaves of frequency 15 GHz are applied. The 
visible lines are the PAT resonances between the left and middle dot near point L in Fig. 3b. 

voltages to energies via a set of conversion factors aj, where i stands for LP, MP or RP 
and j for LQD, MQD or RQD (left, middle and right quantum dot). Here a\ expresses 
by how much a change in the voltage applied to gate i shifts the electrochemical poten- 
tial of dot j. The values of the a^'s can be established using PAT measurements, using 
the known energy of the microwaves as a reference. Figure S3. a shows a measured charge 
stability diagram near the (1, 0, 0)-(0, 1, 0) transition, with microwaves applied. There are 
multiple PAT sidebands visible. The distance between subsequent PAT resonances is set 
by the energy of the microwaves, hv. The slopes of the charging lines of each dot re- 
late the relative influences of each gate on the electrochemical potentials on each dot: 



S LQD = AV RPL /AV LPL = -OL L L Q p D '/ 'a L R Q p D and S M qd = AV RPM /AV LPM ■■ 
Along the charging line of the left dot, the electrochemical potential of the involved state 
on that dot does not change: d/i£,(l,0, 0) = 0. The distance between subsequent PAT reso- 
nances along that charging line then is related to a shift in electrochemical potential of the 
state of the other dot: g?/^m(0, 1,0) = —oi^p D dViPL—otp 1 ®® (IVrpl- Using the expressions re- 



-a 



MQD I MQD 



LP 



jot 



RP 



hv 



-. Similar methods apply for 



lating c^'s to slopes, we can derive: a^% D = . , 
the determination of the other conversion factors. These conversion factors are used for de- 
termining the detuning axis of Fig.3b. The point in gate space where /x.£,(l, 1, 0) = fi M (0, 2, 0) 
serves as a reference point, for which we define 8\ to be zero. When gate voltage i is changed, 
the new value of 8\ is given by the gate voltage change multiplied by a\ — af 1 . A similar 
reasoning applies to 62- 
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VII. FREQUENCY DEPENDENCE OF PAT AND PACT 
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Figure S4: Numerical derivative (along horizontal axis) of the SQD lock-in signal (as in Fig. 4b) 
as a function of frequency versus the detuning between: a (1, 0, 0) and (0, 1, 0), b (0, 0, 1) and 
(0, 1, 0), and c (1, 1, 0) and (0, 1, 1). 

From the frequency dependence of PAT between LQD and MQD (Figure S4.a), and 
between RQD and MQD (Figure S4.b), values for the tunnel couplings between these pairs 
of dots can be established, as in Oosterkamp et al [5]. We find a tunnel coupling strength 
of 12.3//eV between LQD and MQD and of 8.1/xeV between RQD and MQD, by fitting 
using AE = 2a/ [hv) 2 — (2t c ) 2 ). As expected, for high driving frequencies, the frequency of 
the PACT resonances (Figure S4.c) shows a linear dependence on detuning, for each of the 
visible multiphoton resonances. 
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